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We  develop  a  computational  method^ for  the  estimation  of  parame¬ 
ters  in  a  distributed  model  for  a  flexible  structure.  The  structure 
we  consider  (part  of  the  “1*PL  experiment")  consists  of  a  cantilevered 
beam  with  a  thruster  and  linear  accelerometer  at  the  free  end.  The 
thruster  is  fed  by  a  pressurized  hose  whose  horizontal  motion  effects 
the  transverse  vibration  of  the  beam.  We  us^  the  Euler-Bernoulli 
theory f to  model  the  vibration  of  the  beam  and  treat  the  hose- thruster 
assembly  as  a  lumped  or  point  mass-dashpot-spring  system  at  the  tip. 

.  i  f 

-Using,  measurements  of  linear  aooeleratlon  at  the  tip^,  estimate  the 
hose  parameters  (mass,  stiffness,  damping)  and  a  Voigt-Kelvin 
viscoelastic  structural  damping  parameter  for  the  beam  using  a  least 
squares  fit  to  the  data. 

We  consider  spline  based  approximations  to  the  hybrid  (coupled 
ordinary  and  partial  differential  equations)  system;  theoretical 
convergence  results  and  numerical  studies  with  both  simulation  and 
actual  experimental  data  obtained  from  the  structure  are  presented  and 
discussed . 


The  difficulties  involved  in  the  design  of  practical  and  effi¬ 
cient  control  laws  for  large  flexible  spacecraft  (e.g.  the  inherent 
infinite  dimensionality  of  the  system,  a  large  number  of  closely 
spaced  modal  frequencies,  high  flexibility,  light  damping,  a  fuel- 
limited,  hostile,  highly  variable  environment,  etc.)  have  stimulated 
research  into  the  development  of  system  identification  and  parameter 
estimation  procedures  which  will  yield  high  fidelity  models.  A  partic¬ 
ular  area  of  interest  involves  schemes  for  the  estimation  of  material 
parameters  describing,  for  example,  mass,  inertia,  stiffness  or 
damping  properties  in  distributed  models  for  the  vibration  of 
viscoelastic  systems-specifically ,  mechanical  beams,  plates  and  the 
like.  In  addition,  since  the  resulting  inverse  problems  are  often 
infinite  dimensional,  substantial  attention  has  been  focused  on 
approximation;  see,  for  example,  [1],  [2],  [3],  [4],  [ 8 D  and  [12]. 

In  these  treatments,  the  parameter  estimation  problem  is  formulated  as 
a  least  squares  fit  to  measurements  of  either  displacement  or 
velocity.  Although  significant  gains  have  been  made  in  the  development 
of  Instrumentation  to  measure  displacement  and  velocity  (e.g.  laser 
technology,  etc.),  one  of  the  least  expensive,  most  reliable  and  most 
commonly  used  sensors  is  the  linear  accelerometer.  While  in  principle 
it  is  possible  to  integrate  acceleration  measurements  once  or  twice 
to  obtain  respectively  velocity  or  displacement  data,  in  practice  this 
task  oan  pose  significant  challenges.  For  example,  Integration  of  the 
signal  could  result  in  the  amplification  of  low  frequency  measurement 
noise  or  dynamic  effects  which  have  not  been  included  in  the  underly¬ 
ing  model.  In  light  of  this,  we  have  undertaken  to  show  here,  both 


theoretically  and  computationally,  that  a  scheme  in  the  spirit  of 
those  developed  in  the  previously  cited  references  can  also  be 
effectively  used  with  acceleration  measurements.  In  particular  we 
note,  this  involves  the  nontrivial  extension  of  the  familiar 
variational  arguments  which  are  used  to  demonstrate  the  convergence  of 
the  finite  element  state  approximations  upon  which  the  identification 
schemes  are  based.  Indeed,  it  must  be  shown  that  in  addition  to  the 
convergence  of  the  displacement  and  velocity,  the  convergence  of 
acceleration  can  be  obtained  as  well. 

The  other  primary  motivation  for  the  present  effort  is  that  while 
these  methods  have  been  extensively  tested  and  evaluated  with  simula¬ 
tion  data,  they  have  never  been  tried  with  actual  experimental  data. 

We  have  tested  our  scheme  with  data  obtained  from  an  experimental 
structure  which  was  designed  and  constructed  at  the  Charles  Stark 
Draper  Laboratory  in  Cambridge,  Massachusetts  with  funding  provided  by 
the  United  States  Air  Force  Rocket  Propulsion  Laboratory  (RPL) .  The 
RPL  structure  (as  it  will  henceforth  be  referred  to  as)  was  designed 
to  serve  as  a  test  bed  for  the  implementation  and  evaluation  of 
control  algorithms  for  large  angle  slewing  of  spacecraft  with  flexible 
appendages.  The  structure  was  specifically  designed  to  exhibit 
structural  modes  and  damping  characteristics  representative  of 
realistic  large  flexible  space  structures. 

In  Section  2  we  describe  the  RPL  structure  (its  geometry, 
instrumentation,  etc.)  and  formulate  an  inverse  problem  involving  a 
distributed  system.  In  Section  3,  we  use  the  resulting  infinite 
dimensional  estimation  problem  to  motivate  the  development  of  a  finite 
dimensional ,  finite  element  based  approximation  scheme .  We  also 


discuss  our  theoretical  convergence  results.  In  Section  4  we  present 
numerical  findings. 

We  use  standard  notation  throughout.  For  X  a  normed  linear 
space,  L(X)  denotes  the  space  of  bounded  linear  operators  from  X  into 

lr 

X.  For  ft  an  interval  and  k  -  0,1, 2,*-*,  C  (ft;X)  denotes  the  space  of 
functions  from  ft  into  X  which  are  k  times  continuously  strongly 
differentiable  on  ft.  When  k  -  0  we  shall  simply  write  C(ft;X).  A 
function  f  from  ft  into  X  will  be  said  to  belong  to  L2(ft;X)  if 
Jft|f(t)|^  dt  <  oo.  For  k  -  0, 1 ,2,  •  •  • ,  Hk(ft;X)  denotes  the  completion 

V 

of  C  (ft;X)  with  respect  to  the  norm 

M*  -  (J0  In  |fC3)Ct)|;><n)1/2. 

If,  in  addition,  X  is  a  Hilbert  space  with  inner  product  <  •  ,  •  >x  then 

Jr 

H  (ft;X)is  a  Hilbert  space  with  inner  product 

<f.g>k  “  lQ  g(J)(t)>x  dt. 

When  X  -  R,  we  use  the  abbreviated  notations  Ck(ft),  L,(ft)  and  Hk(ft). 
Note  that  H#(ft)  -  L„(ft)  and  is  the  standard  inner  product  on 


overall  symmetry  of  the  structure.  Nitrogen  gas  from  tanks  mounted  on 
the  central  hub  is  supplied  to  the  thrusters  via  two  stainless  steel 
mesh-wrapped  high  pressure  hoses.  The  expulsion  of  propellant  from 
the  thruster  nozzles  is  controlled  by  electro-mechanical  or  solenoidal 
valves.  Each  of  the  four  appendages  is  equipped  with  a  sensor  in  the 
form  of  a  linear  accelerometer  attached  at  its  tip.  Data  from  the 
accelerometers  is  processed  and  recorded  and  control  input  signals  to 
the  thrusters  are  generated  by  a  MINC  11/23  microcomputer.  A  detailed 
description  of  the  structure's  design  specifications  can  be  found  in 
[6]  and  [ 15] . 

The  problem  which  is  of  primary  concern  to  us  here  involves  the 
modeling  of  the  effects  of  the  nitrogen  supply  hoses  on  the  transverse 
vibration  of  the  active  members.  We  consider  therefore,  the  structure 
with  the  central  hub  Immobilized  and  look  only  at  the  vibration  of  one 
of  the  active  appendages  and  view  it  as  a  simple  cantilevered  beam 
(see  Figure  2.2). 


Figure  2.2 


we  treat  the  thruster  assembly  as  a  point  mass  that  is  rigidly  attach 
ed  to  the  beam  at  the  tip  and  propose  a  model  for  the  hose  effects  in 
the  form  of  a  proof  mass  which  reacts  against  the  tip  mass.  In 
effect,  we  consider  the  idealized,  simplified  structure  depicted  in 
Figure  2.3  below  involving  a  single,  cantilevered,  flexible,  uniform 
beam  with  a  two-mass-dashpot-sprlng  system  affixed  to  its  free  end. 


|ya«iaa^l 


IHKODCIK^I 


Figure  2.3 


In  formulating  a  mathematical  model  for  the  structure  shown  in 
Figure  2.3  above,  we  assume  that  the  beam  is  of  length  a,  with  uniform 
rectangular  cross  section  of  height  h  and  width  b.  We  let  u(t,x)  and 
y(t)  denote  respectively  the  transverse  displacement  of  the  beam  at 
position  x  along  its  span  and  the  displacement  of  the  proof  or  hose 
mass,  each  at  time  t.  Both  are  measured  relative  to  the  x-axis  in  the 
coordinate  frame  determined  by  the  longitudinal  axis  of  the  beam  in 


its  unde formed  state  with  origin  located  at  the  beam's  root  or  fixed 
end.  Assuming  the  beam  undergoes  only  small  deformations  (i.e. 

3m 

|u(t,x)|  <<  X.  and  |— (t,x)|  <<  1)  and  has  a  small  height  to  span 

length  ratio,  the  Euler-Bernoulli  theory  (see  [5])  including  Voigt- 
Kelvin  viscoelastic  structural  damping  (see  [10])  yields  the  partial 
differential  equation 


9*u  3*  3u  94u 

(2.1)  p  - (t.x)  +  c^I - (t,x)  +  El  - (t,x) 


9t; 


9x4  St 


Sx4 


=  0, 

0  <  x  <  8. , 


t  >  0 


where  p  is  the  linear  mass  density  of  the  beam,  E  is  the  modulus  of 
elasticity,  cD  is  the  coefficient  of  viscosity  and  I  is  the  second 
moment  or  moment  of  inertia  of  the  cross  sectional  area  A  about  the 
neutral  axis.  For  the  beam  we  consider  here  with  constant  rectangular 
cross  section,  I  =  bh  /12.  Since  the  beam  is  assumed  to  be  uniform, 
the  parameters  p,  E  and  cQ  are  taken  to  be  constant  in  time  and  space. 

Balancing  forces  at  the  free  end,  elementary  Newtonian  mechanics 
yields  the  equations  of  motion 

9*u  S3  9u,  33u 

(2.2)  m_  - (t , it)  -  cnl - (t,i)  -  El  - (t,X.) 

1  9t 1  9x3  3t  3x3 

dv  9u 

-  c„  (— (t) - (t,Jl))  +  k„  (y(t)  -  u(t,D)  +  f(t),  t  >  0 

H  dt  9t  H 

and 

d*v  dv  3u 

(2.3)  m„  — -(t)  +  c„  (— (t) - (t.X.))  +  k„  (y(t)  -  u( t , X. ) )  -  0, 

H  dt*  H  dt  3t  H 


t  >  0 


for  the  tip  and  hose  masses  mT  and  mH  respectively.  Here  kH  is  the 
hose  stiffness,  cH  is  the  hose  damping  coefficient  and  f(t)  is  the 
externally  applied  force  at  time  t  due  to  the  firing  of  the  thrusters 
mounted  at  the  tip. 

Making  the  assumption  that  the  rotatory  inertia  of  the  proof  mass 
system  is  negligible,  rotational  equilibrium  at  the  tip  can  be 
expressed  as 

91  9u  8*u 

(2.4)  cnI - (t,!l)  +  El  - (t.JL)  =  0,  t  >  0. 

u  ax1  at  ax* 

The  zero  displacement  and  zero  slope  constraints  at  the  fixed  end  are 
given  by 

3u 

(2.5)  u(t,0)  -  0  and  — (t,0)  -  0,  t  >  0 

8x 

respectively.  Taking  the  structure  to  be  initially  at  rest  we  have 
the  Initial  conditions 

(2.6)  u(0,x)  «=  0  and 

and 

(2.7)  y(0)  -  0  and 


au 

— (0,x)  =  o, 

at 


0  £  x  £  51 


dy 

i(0)  - 


In  the  mathematical  model  given  by  (2.1)  -  (2.7)  above  the  parameters 
p,  mT  and  I  can  be  measured  or  oomputed  directly.  The  modulus  of 
elasticity  E  is  typically  determined  in  the  laboratory.  For  the  most 
commonly  used  materials  (including  aluminum  which  is  the  material  from 


which  the  structure  of  interest  to  us  here  is  made)  its  value  can  be 
readily  looked  up  in  standard  engineering  tables.  The  parameters  cD, 
mH*  °h  and  kH  on  the  otlier  hand,  must  be  determined  experimentally; 
that  is,  they  will  have  to  be  identified  based  upon  the  observed 
response  of  the  structure  to  a  given  input  disturbance.  This  is  one 
class  of  inverse  problems  which  we  formulate  and  consider  below.  In 
the  system  of  equations  (2.1)  -  (2.7)  we  explicitly  modeled  (albeit, 
in  a  rather  simple  fashion)  the  dynamical  effects  of  the  hose.  The 
unknown  hose  parameters  are  then  determined  as  the  solution  to  an 
inverse  problem. 

An  alternative  approach  to  obtaining  a  model  which  exhibits  a 
reasonable  degree  of  fidelity  involves  a  technique  which  is  sometimes 
referred  to  as  model  adjustment.  Starting  with  a  simple  model,  the 
parameters  are  then  "adjusted"  so  as  to  compensate  for  unmodeled 
dynamics.  The  choice  of  parameters  to  be  adjusted  and  the  resulting 
variations  may  or  may  not  be  motivated  by  physical  considerations. 

In  our  problem  for  example,  we  might  consider  a  simple  cantile¬ 
vered  beam  with  tip  mass  (i.e.  mH  -  cH  -  kR  -  0)  and  then  adjust  the 
theoretical  or  measured  values  of  E  and  mT  to  compensate  for  the 
dynamical  effects  which  result  from  the  hose  mass  and  motion.  A  value 
for  the  parameter  cD  could  also  be  identified  if  damping  effects  are 
considered  significant.  Model  adjustment  was  used  in  [6]  to  obtain  a 
model  for  the  RPL  structure  upon  which  control  design  could  be  based. 

We  define  an  inverse  problem  which  encompasses  both  of  the 
general  approaches  which  have  been  outlined  above.  We  assume  that  an 
input  disturbance  described  by  the  function  f(t),  t  e  [0,T]  is  applied 
to  the  structure  via  the  tip  thrusters  and  that  the  linear  accelera- 
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tion  at  the  free  end  of  the  beam,  z(t),  is  measured  and  recorded  for 
each  t  €  [to.ti]  where  0  £  ta  £  £  T .  (Of  course,  in  actual 

practice,  z  could  in  fact  only  be  sampled  discretely).  Let  R+  denote 
the  positive  real  numbers  and  let  Q  be  a  closed  and  bounded  subset  of 
R°.  We  seek  a  q  €  <}  which  minimizes 

JCq)  -  f^1  |^-^(t,iL;q)  -  z(t)|2dt 
1 1 1  1 3t 1  1 

where  u(-,*;q)  denotes  the  solution  to  the  initial-boundary  value 
problem  (2.1)  -  (2.7)  corresponding  to  q  -  (m^, , E , c^ , m^ , cH , kH )  €  Q. 

Our  primary  concerns  in  the  next  section  will  include  well- 


posedness  of  the  system  (2.1)  -  (2.7),  existence  of  a  mlnimizer  for  J, 
and  development  of  approximation  techniques  to  find  this  minimizer. 


A  computational  method  for  the  solution  of  the  estimation  problem 
posed  above  will  invariably  involve  finite  dimensional  approximation 
of  the  initial -boundary  value  problem  (2.1)  -  (2.7).  We  have  been 
successful  in  solving  Inverse  problems  for  distributed  parameter 
models  for  flexible  structures  (see,  for  example,  [1],  [2],  [3],  [4], 
t 12] )  using  spline-based  Ritz-Galerkin  techniques.  We  apply  those 
ideas  here  and  derive  finite  element  approximations  based  upon  an 
abstract  Hilbert  space  formulation  of  the  hybrid  system  of  ordinary 
and  partial  differential  equations  and  boundary  conditions  given  in 
(2.1)  -  (2.7).  This  abstract  formulation  is  also  useful  in 
establishing  existence,  uniqueness  and  necessary  regularity  results 
for  solutions.  We  briefly  outline  the  essential  features  of  our 
general  approach  (Including  theoretical  convergence  results)  in  the 
context  of  the  particular  problem  of  interest  to  us  here. 

Let  H  -  R*  x  L, (0 , i)  be  endowed  with  the  usual  product  space 
inner  product 

<($.ti,$).(x,u,nO>h  -  $x  +  mi  +  < > 0 

and  let 

V  -  {($,71,$)  €  H:  $  €  H*(0,JO,  $(0)  -  D$(0)  -  0,  71  -  $(*.)} 

be  endowed  with  the  inner  product 

<($,$(*.),$),  (X  .*(*.), v)>v  -  ($-$(S.))(X-N»0O)  +  <D*$,D,V>0 


12 


where  the  symbol  D  is  used  here  and  below  to  denote  the  spatial 

d 

differentiation  operator  — .  The  space  V  together  with  the  inner 

product  form  a  Hilbert  space  which  is  densely  and  compactly 

embedded  in  H. 

We  rewrite  the  system  (2.1)  -  (2.7)  as  the  abstract  second  order 
initial  value  problem  in  H 

(3.1)  Mutt(t)  +  Cut(t)  +  Ku(t)  -  F(t),  t  >  0 

(3.2)  S(u(t)  +  eut(t))  -  0,  t  >  0 

(3.3)  u(0)  -  0  ut(0)  -  0 

in  the  states  u(t)  -  (y(t) ,u(t,  JQ  ,u(t, • )) .  The  operators  M  €  L(H) , 
C:D  c  H  -♦  H  and  K:D  c  H  -»  H  are  given  by 

(3.4)  M(£  ,Ti»<t>)  -  (mH^  .m^ri,  pO)  , 

C(S.t,40  -  (cH(^-u),cH(t-U  -  oDl  D*<|)(!0,cdI  D«<t>), 

and 

K(?,ti,<1>)  -  (kH(^-Tl),kH(n-U  -  El  D*0(Jt)  , El  D * 0 ) 

where  D  -  ,ti,0)  e  V:  <t>  €  H4(0,i)>.  For  each  t  >  0,  F(t)  - 

(0,f(t),0)  €  H,  6:  D  c  H  -♦  R  is  given  by  6(U,n,40)  -  D*4>(0  and  e  - 
cd/E. 

The  restrictions  C  and  K  of  the  operators  C  and  K  that  appear  in 
equation  (3.1)  above  to  N(6),  the  null  space  of  the  operator  6,  have 


7? 


natural  extensions  to  bounded  operators  from  V  (which  is  the  V-closure 
of  N(5))  into  V',  the  dual  of  V.  The  extensions  are  defined  in  terms 
of  the  bilinear  forms  c(  •  ,  •  ) :  V  x  V  -»  R  and  k(  •  ,  • ) :  V  x  V  -»  R  given  by 


(3.5)  (C<fr)(iO  -  c(4>,n>)  -  ch(;-4>(JO)(X-n»0O)  +  cDI<D‘4),DlS'>0 
and 

(3.6)  (kJ)(v)  -  k(i,i)  -  kH(?-^(0)(X-lr(JL))  +  EI<D*<|>,DlV>0 

for  i  -  (2; , d>(«.) , 0)  €  V  and  v  -  (X. ♦(!),♦)  e  V. 

The  finite  element  method  we  develop  below  could  be  derived  from 
standard  energy  considerations .  While  this  is  not  the  approach  we 
take,  it  is  worth  noting  that  the  usual  energy  expressions  can  be 
given  in  terms  of  the  forms,  operators  and  inner  products  defined 
above.  The  kinetic  energy  is  given  by 

1 

T0  ’  T  <Mut(t),ut(t)>H, 

the  potential  or  strain  energy  by 
1 

Un  -  —  k(u(t) ,u(t) ) 
u  2 

and  the  Rayleigh  dissipation  function  by 

Fq  -  —  c(ut(t) ,ut(t)) . 

2 

Written  in  its  weak,  variational  or  distributional  form 
<Mu..  (t)  ,<J>>W  +  c(u.  (t)  ,<}>)  +  k(u(t)  ,  $)  -  <F(t)  .  <t>>TT, 


(3.7) 
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(3.8)  u(0)  -  0  ut(0)  “  0 

the  initial  value  problem  (3.1)  -  (3.2)  in  H  becomes  an  initial  value 
problem  in  V'.  If  we  assume  that  f  €  L,(0,T)  and  rewrite  (3.7).  (3.8) 
as  an  equivalent  first  order  vector  system,  the  theory  of  abstract 
parabolic  systems  (see  [93,  [143)  yields  the  existence  of  a  unique 
mapping 

u  6  C([0,T] ;V)  n  H1 ( (0 , T) ; V)  H  Cx([0,T];H)  O  H1((0,T);V') 

which  satisfies  (3.7),  (3.8).  If  we  are  willing  to  assume  further 
that  f  is  Holder  continuous  then  there  exists  a 

(3.9)  u  €  C([0,T] ;V)  n  CX((0,T3;V)  n  CX([0,T3;H)  n  C*((0,T];H) 

A  A 

with  u(t)  +  eut(t)  6  D,  t  >  0  which  uniquely  satisfies  the  initial 
value  problem  (3.1)  -  (3.3). 

In  order  to  demonstrate  the  convergence  of  the  approximation 
schemes  we  develop  below,  we  shall  require  a  somewhat  more  regular 
solution  to  the  initial  value  problem  (3.7),  (3.8)  than  either  of  the 
conditions  on  f  stated  above  can  guarantee.  In  addition  to  (3.9),  we 

a 

shall  require  that  u  €  H*((0,T);V).  This  can  be  guaranteed  (see  [7]) 
if  we  assume  that  f  €  Hx(-t,T)  for  some  t  >  0  with  f(t)  -  0,  t  <  0  and 
we  modify  our  original  mathematical  model  so  that 

(3.10)  F(t)  -  f(t)0,  t  €  [-t,T] 

for  some  8  -  (0,0(1), 0),  a  fixed  element  in  V.  We  note  that  with  0 


1 


f, 
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chosen  appropriately  in  V,  F  given  by  (3.10)  may  in  fact  represent  an 
improved  model  of  reality  when  compared  with  our  present  choice  of  F 

a 

where  0  -  (0,1,0)  e  H. 

Central  to  our  approach  is  a  cubic  spline  based  Galerkin  approxi¬ 
mation  to  the  initial  value  problem  (3.7),  (3.8).  For  each  N  «1,2,--- 


N  l  2JI 

let  A  denote  the  uniform  mesh  {0,  — ,  — ,  •••,  Jl}  on  [0,51]  and  let 

N  N 

r  N  N+l 

{B  >  denote  the  usual  cubic  B-splines  defined  with  respect  to  the 

3  3—1 

nodal  set  A  (see  [11],  [13]).  Briefly,  each  BN  is  a  function  on 

3 

jt  g. 

1 0 ,  JL ]  which  is  a  cubic  polynomial  on  each  subinterval  [  (k-l>^-,k—  ] , 

N  r  l  l, 

k  -  1,2,  •  •  •  ,N .  The  support  of  B^  is  ^  C  J-2)^-,  (J+2)^-]  n  [ 0 ,  JL  ]  with 

B*(£)  “  4*  DB/3^)  “  0,  b*((3±i£)  -  l.  and  db*((3±1)4)  -  + 

3  N  JN  jN  J  N  JL 

Defining  (6  }  by  e  -  B  -2B  -2B_  and  3  -  B  ,  3-2,3, ••• ,N+1 , 

J  w  ^  10  1  1  J  J 

N  x  N 

we  have  3  (0)  -  DB  (0)  -  0,  j  -  1,2, • • • .N+l.  With  3  -  (1,0,0)  and 

3  J  0 

N  N  N  N  N+l 

3  -  (0, 3  .(Jl) ,  3  ) ,  3  “  1,2,***,N+1,  V  -  span  (3  }  is  an  N+2 

333  33-0 

dimensional  subspace  of  V. 

The  Galerkln  equations  in  VN  corresponding  to  (3.7),  (3.8)  for 
~N  %  N 

u  (t)  €  V  are  given  by 


-N  aN  .N  aN  aN  .N  .N 

(3.11)  <Hu  (t ) , 3  >  +  c(u  (t ) , 3  )  +  k(u  (t ) , 3  )  -  <F(t),8  >  , 

tt  J  H  t  j  j  JH 

t  >  0,  3  -  0,1,2, • . • , N+l 

~N  aN 

(3.12)  u  (0)  -  0  u  (0)  -  0. 

ti 


Setting 

N+l  N  ~N 

u  (t)  -  T,  w  (t)6  ,  t  *  0, 

J-0  J  J 

N 

the  initial  value  problem  (3.11),  (3.12)  in  V  is  equivalent  to  the 
linear,  nonhomogeneous ,  second  order  N+2  -  vector  system 


(3.13)  M 


N 


d*w 

dt' 


dw 


,N 


-(t)  +  CN  ^-(t)  +  KNwN(t) 
dt 


(3.14)  wN(0)  -  0 


dw 


,N 


dt 


-(0) 


N  N  N  NT 

where  w  (t)  -  (w  (t),w  (t),---,w  (t))  . 

0  1  N+l 

.  .  N  N  N 

x  (N+2)  matrices  M  ,  C  and  K  are  given  by 


and 


N 

& 

< 

< 

M 

-  <M0  ,6  > 

i.J 

i  J  H 

N 

Jz; 

< 

& 

< 

C 

-  c( 8  ,8  ), 

i.J 

i  J 

N 

~N  ~N 

Ki.j ' 


FN(t),  t  >  0 


0 


The  entries  in  the  (N+2) 


i.J  “  0,1,2, ••• ,N+1  respectively.  For  each  t  >  0  the  components  in 

N  N  „N  N 

the  N+2  -  vector  F  (t)  are  given  by  F  (t)  -  <F(t),8  >  -  f(t)8  (i) 

i  i  H  i 


or,  recalling  (3.10),  by 


p  (t)  -  f (t) <0,6  >  -  f(t)(eOL)e  co  +  <e,e  >  ), 

i  i  H  i  i  0 


i  -  0,1,2, •• • ,N+1 . 


We  consider  the  sequence  of  approximating  finite  dimensional  iden¬ 
tification  problems  which  consist  of  finding  qN  6  Q  which  minimizes 


(3.15) 


.N 


Aq)  -  J 


2 

dt 


where  for  each  q  £  Q,  uN(t;q)  -  (yN(t;q)  ,uN(t .  JL  ;  q)  ,uN(t ,  •  ;q))  is  the 
unique  solution  to  the  initial  value  problem  (3.11),  (3.12)  in  VN 
corresponding  to  q  -  (mT,E,cD,mH>cH,kH)  e  Q.  In  actual  practice,  for 
a  given  q  €  Q,  JN(q)  is  computed  as 


JN(q)  -  lWN-i(t:q)  +  4wN(t;<*> 


N  , 

+  w  (t;q)  -  z(t) 
N+l 


dt 


N  N  NT 

where  w  ( - ; q)  -  (w  ( - ; q) , *  *  - , w  (*;q))  is  the  unique  solution  to  the 

0  N+l 

N+2  -  vector  system  (3.13),  (3.14)  corresponding  to  q  €  Q. 

With  finite  dimensional  state  constraints,  the  solution  of  the 
N  estimation  problem  above  is,  at  least  in  principle,  routine.  For 
inverse  problems  which  are  closely  related  to  the  one  we  treat  here, 
our  earlier  numerical  studies  have  shown  that  satisfactory  results  can 
be  obtained  using  any  one  of  a  number  of  standard  computational 
techniques  for  least  squares  minimization  (for  example,  Newton's 
method,  conjugate  gradient,  steepest  descent,  Levenberg-Marquardt , 
etc. ,  see  [2] ) . 
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Our  fundamental  theoretical  result  is  that  each  of  the  approxi¬ 
mating  identification  problems  and  the  original  problem  have 
solutions.  Moreover,  we  show  that  the  solutions  to  the  approximating 
problems,  in  some  sense,  approximate  solutions  to  the  original 
problem.  We  require  the  following  lemma. 

T-emma.  s.i  Suppose  {qN}  C  Q  with  qN  -»  q°  as  N  -*  Let  uN(  •  ;qN) 

denote  the  unique  solution  to  the  initial  value  problem  (3.11),  (3.12) 

N  ~  n 

corresponding  to  q  and  let  u(-;q  )  denote  the  unique  solution  to  the 
Initial  value  problem  (3.7),  (3.8)  corresponding  to  q°.  If  u(-;q°)  € 
H1 ( (0 , T) ; V)  then 


as  N  -»  oo. 

Proof 

For  each  N  «  1,2, •••  let  PN  denote  the  orthogonal  projection  of  H 
onto  VN  defined  with  respect  to  the  standard  inner  product  on  H, 

Using  the  approximation  theoretic  properties  of  interpolatory 
splines,  it  Is  not  difficult  to  show  that  (see  [3]) 

(3.17)  11m  |(PN  -  I)U, 71,0)1  -  0 

N-*oo  i  •  H 

for  each  (^,^,0)  €  H  and  that 

(3.18)  lim  |(PN  -  I)<T> I  - 


0 


for  each  $  €  V. 


For  q  -  (mT>E,cD,mH,cH,kH)  e  Q  it  is  immediately  clear  that  M , 
c(-,-)  and  k( • , • ) ,  the  operator  and  forms  defined  in  (3.4),  (3.5)  and 

(3.6)  respectively  depend  upon  q.  For  q°  =  (m0,E0,c0,m0,c0,k0)  e  Q 

T  D  H  H  H 

N  ✓NNNNNN 

and  q  -(m,E,c,m,c,k)€Qwe  adopt  the  shorthand  notation 
T  DHHH 

M°  -  H(q°).  C°( • , • )  -  c(q°)( • , • ) ,  k°(-.-)  -  k(q°)(-,-),  MN  -  M(qN) , 
cN( • , • )  -  c(qN)( • , • )  and  kN( • , • )  -  k(qN)( • , • ) .  Similarly,  we  denote 

r\  ^  >V  -lyT  * 

u(*;q  )  and  u  (*;q  )  by  u  and  u  respectively. 

From  (3.17),  the  assumption  that  u°  6  H*((0,T);V)  and  the 
inequality 


it  is  clear  that  we  need  only  to  consider  the  first  term  on  the  right 
hand  side  of  the  above  estimate . 

Letting  vN(t)  -  uN(t)  -  PNu°(t)  for  t  s  0,  (3.7),  (3.8),  (3.11), 
(3.12)  and  VN  C  V  imply 


(3.19) 


{Mn;n  :n  „n,,;n  :n.  .  .n^n  :n. 

vtt’^  H  +  °  ^Vt’^  )  +  k  (V  ,<t»  ) 


-  ♦  <(M°  - 

*  oN(«-pH)u°,i1')  ♦  -  0N(G°,iH) 


.h,,,.  „n,-o  :n,  .  ^o,:o  :n,  ,_n,-o  :n 


+  k"((I-P")uu,0  +  kuCu",0  -  k"(uu,*").  t  >  o,  €  vM 


Choosing  4>  -  v^Ct)  G  V  ,  from  (3.19)  we  obtain 

“  <HNCl-PN)u?t.v;t>H  +  <CH°-HN)^t.v?t>H 

+  ft  CN((I-PN)u°,^)  -  CN((I-PN)^t.^) 

+  St  ~  '  <O0(u°t.^)  -  0N(^t.v?)} 

+  d_kN((I_pN);0>;N)  _  kN((I_pN)G0  ^ 

+  g_  {kO(;Ot;N)  .  kN(;o  ;N)}  _  {ltO£0'fy  _  kN(;o  ;N)} 


-ft  k*<^>  ♦  kNc^,;N)( 


t  >  c, 


Integrating  the  above  expression  from  0  to  t  and  recalling  (3.20),  we 


find 


(3.21) 


-  J0  <’«N«-^E°S-LV  ■«0-hn^°s-L-h 

-  oNca-pi,)i°s.v»)  -  (o°c;°s.^)  -  oN(;°s.^)) 

-  *KCCI -PH)i°.^)  -  (k°(i°.v“)  -  kN(i°.v^) 

+  *HC^*))dS 

♦  cnc«-pn);°  o?)  ♦  (o°c;°.v?)  -  oN(;?.i?)) 
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+  kN((i-Pi,)G°,v»)  ♦  (k°ci°,v^  -  kK(;°,^))  -  kNcvN,^). 

We  recall  that  Q  has  been  assumed,  to  be  a  closed  and  bounded  subset  of 
R®  and  observe  therefore  that  the  forms  c°(-,-),  cN(-,-).  k°( • , • )  and 
k  (•,•)  are  uniformly  bounded.  These  two  facts  together  with  the 
repeated  application  of  the  inequality 


<a,b>  S  | a |  1  b |  £  a|a|2  +  -  | Id | 2 ,  a  >  0 

4a 


in  (3.21)  yield  the  estimate 


i  ivLi: 


ds  +  V 


s  v  i  &  ici-^2.1!  *  ■rc.i: 


I  1  r  lrnN  rr,0!2  ■  |JLJ>|8 

4a  M  H  mH I  I “t  roT I 


I  +  h?-”Sl  H^ssl!  ♦ 


♦  AIL ♦  ♦  Ms^Ll!  ♦  K\ 


+  *  Wi:  -  <M" ♦  I^V^O  tfl 


*  *  |c^i2>|u?C  *  °ri;  * 


♦  -iv?i2  ♦  i-cik2-k°i2  ♦  ien-e0i2)i;°i2  +  a i i 2 
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+ 


+ 


where  Yq  is  a  positive  constant. 


Choosing  a  >  0  sufficiently  small , 


we  find 


where 


and  Yi,  i  -  1,2,3  are  positive  constants  which  do  not  depend  on  N. 

Alt  A  IT  M 

Choosing  0  -  v^(t)  6  V  in  (3.19),  arguments  similar  to  those 

used  above  (see  [2],  [3])  yield 

(3.23)  lim  |vN(t)|2  -  0 

N-®  1  1 v 

for  each  t  e  [0,T].  Using  u°  €  H*((0,T);V),  (3.18)  and  an  application 
of  the  Gronwall  Inequality  to  (3.22)  we  obtain  the  desired  result. 


We  note  that  we  also  obtain 


.  A  . ^  *  ■  vut. r. 


(3.24) 


lim  v"(t) 
N^®  1  z  !V 


for  each  t  e  [0,T].  From  (3.23)  and  (3.24)  we  find  uN(t ;qN)-u(t ;q°) 


-*  0  and  |u^(t;qN)-ut(t  ;q°)  |  -»  0  as  N  -*  »  for  each  t  €  [0,T]. 

We  remark  that  it  is  the  L»  convergence  (more  precisely,  H 
convergence)  in  (3.16)  which  necessitates,  at  least  in  theory,  that  we 
be  provided  with  distributed  time  observations  (i.e.  observations 
which  are  continuous  in  time).  It  is  clear  from  (3.23)  and  (3.24) 
that  for  fits  based  upon  displacement,  velocity  or  slope,  time-sampled 
measurements  are  sufficient.  Of  course  when  the  approximating 
optimization  problems  are  solved,  the  integral  least  squares 
performance  indices  (3.15)  are  discretized.  Consequently,  in 
practice,  only  discrete  measurements  of  linear  acceleration  at  the  tip 
are  required. 


Theorem  3.1  Each  of  the  approximating  identification  problems  has  a 
solution  q  .  The  sequence  {q  }  c  Q  admits  a  convergent  subsequence 
-N1  -N1  - 

{q  J}  with  q  J  -*  q  €  Q  as  J  ■*  ».  If  for  each  q  e  <},  u(  •  ;q),  the  unique 
solution  to  the  initial  value  problem  (3.7),  (3.8)  corresponding  to  q, 
is  an  element  in  H*((0,T);V)  then  q  is  a  solution  to  the  original 
identification  problem.  In  addition,  the  limit  point  of  any  convergent 
subsequence  of  {q  }  is  a  solution  to  the  original  identification 
problem  as  well. 

Proof  Standard  continuous  dependence  results  for  linear  ordinary 
differential  equations,  the  fact  that  Q  has  been  assumed  to  be  a 
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dosed  and  bounded  subset  of  R®  and  the  form  of  are  sufficient  to 

oonolude  that  a  solution  qN  €  Q  to  the  Nth  approximating  identifica¬ 
tion  problem  exists.  Once  again  since  Q  is  a  closed  and  bounded  (and 

therefore  compact)  subset  of  R  ,  the  sequence  {q  }  c  Q  admits  a 

_N .  _N  _N, 

convergent  subsequence .  If  (q  C  (q  }  with  qJ-*q€QasJ-*®  and 
q  is  any  point  in  Q,  then  two  applications  of  Lemma  3.1  (the  second 
one  with  the  constant  sequence  {q})  yield 

N,  _N j  H. 

J(q)  -  11m  J  J(q  J)  S  lim  J  J(q)  -  J(q) 

j-*®  j-*® 

and  the  theorem  is  proved. 

Although  Theorem  3.1  above  guarantees  only  subsequential 
convergence,  in  all  test  and  simulation  examples  we  have  considered, 
we  in  fact  observe  the  convergence  of  the  sequence  (qN)  itself  to  the 
optimal  parameters  q.  Also,  it  is  not  difficult  to  verify  that  with 
only  minor  modification  (see  [2])  the  approximation  scheme  reported  on 
here  (together  with  the  convergence  theory  outlined  in  the  lemma  and 
theorem  above)  could  be  applied  to  inverse  problems  involving  the 
estmatlon  of  spatially  varying  parameters  (such  as  linear  mass  density 
p,  flexural  stiffness  El,  or  damping  coefficient  oDI)  which  appear  in 
the  equations  (2.1)  -  (2.4).  We  note  of  oourse  that  when  either  El  or 
cDI  are  spatially  varying,  the  Euler-Bernoulli  equation  and 
corresponding  boundary  conditions  are  of  a  slightly  different  form 
than  those  given  in  (2.1)  -  (2.4)  (see  [3]). 


We  used  our  scheme  to  attempt  to  solve  the  inverse  problem  which 
was  posed  above  with  data  obtained  from  an  experiment  on  the  RPL 
structure.  We  report  on  our  findings  and  observations  here. 

All  computer  codes  were  written  in  Fortran  and  run  on  the  IBM 
3081  at  the  University  of  Southern  California.  The  approximating 
finite  dimensional  least-squares  minimization  problems  were  solved 
using  the  IMSL  implementation  of  the  Levenberg-Harquardt  algorithm 
(routine  ZXSSQ),  an  Iterative  Newton's  method-steepest  descent  hybrid 
(see[2l).  The  second  order  N+2  -  vector  systems  (3.13),  (3.14)  were 
solved  (integrated)  in  each  iteration  (for  the  evaluation  of  JN  and 
its  gradient)  using  Gear's  method  for  stiff  systems  (IMSL  routine 
DGEAR) .  The  integral  least  squares  performance  index  was  approximated 
by  a  discrete  sum  over  a  uniform  mesh  on  It,,tx].  The  integral  inner 
products  in  the  definitions  of  the  matrices  MN,  CN  and  KN  were 
computed  using  a  composite  two  point  Gauss-Legendre  quadrature  rule. 
The  second  time  derivative  of  wN,  or  generalized  acceleration, 

d*wN 

— — ,  was  computed  using  a  second  order  centered  difference  on  the 
dt* 

generalized  displacement , 


(4.1) 


wN(t+A)  -  2wN(t)  +  wN(t-A) 


We  found  this  to  be  a  somewhat  more  stable  method  for  computing 
acceleration  (an  unbounded  measurement)  than  was  a  first  order 
centered  difference  on  the  generalized  velocity. 
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(4.2) 


dvN,  AN  dwN,  A 

- (t  +  -) - (t  -  -) 

dt  2  dt  2 


Either  of  the  time  differencing  formulas  (4.1)  or  (4.2)  proved  to  be 
significantly  more  stable  than  using  the  differential  equation  (3.13) 

d'WH 

directly  to  compute  — (t)  via  an  inversion  of  M  .  As  to  why  this 

was  so,  we  can  only  offer  the  conjeoture  that  the  time  differencing 
provided,  at  least  to  a  certain  extent,  some  filtration  of  the  signal. 

Before  turning  our  attention  to  the  experimental  data,  we  tested 
our  scheme  with  simulated  data.  "True"  values  for  the  unknown 
parameters  cD  (actually  oDI),  mH,  cH  and  kg  were  chosen  and  a  qulntic 
spline-based  semi -discrete  Galerkin  scheme  applied  to  the  initial 
value  problem  (3.7),  (3.8)  was  used  to  generate  data. 

Setting  p  .03,  m  .15,  El  -  80.0,  i  -  4.0  and 


f(t)  -  { 


0  i  t  S  0.05 
0.05  <  t  S  5.0, 


the  fit  was  carried  out  based  upon  observations  of  linear  acceleration 
at  the  tip  at  times  ^  -  .11,  1  -  2, 3, •••,50.  We  note  that  this  is 
equivalent  to  taking  t,  -  .1,  ti  -  5.0  and  using  a  standard  rectangle 
rule  with  uniform  mesh  spacing  . 1  to  discretize  the  integral  appearing 
in  the  definition  of  the  least  squares  performance  index  J1* .  The 
Initial  estimates  oDI  -  .0035,  mg  -  .035,  and  kg  -  .4  were  used  to 
start  the  iterative  optimization  procedure.  In  (4.1),  A  was  taken  to 
be  .1.  Our  results  are  summarized  in  Table  4.1  below. 


rj 


2? 


N 

cNI 

CDI 

-N 

CN 

CH 

kN 

A<f) 

2 

.037537 

.039471 

.003428 

. 298626 

2.57xl0_1 

3 

. 066997 

.039485 

.003907 

.298875 

4.37xl0“2 

4 

.005063 

.039777 

. 003997 

.299455 

5.06xl0“3 

5 

.005667 

.039899 

.003971 

.299787 

7.66xl0~4 

6 

.005049 

.040035 

. 004006 

.300087 

4.63xl0~5 

True 

value 

.005000 

.040000 

.004000 

.300000 

Initial 

Estimate 

.003500 

.035000 

.003500 

.400000 

Table  4 . 1 


The  experiment  which  we  describe  below  was  carried  out  for  us  on 
the  RPL  structure  by  Dr.  Michel  A.  Floyd,  formerly  of  the  Control  and 
Flight  Dynamics  Division  of  the  Charles  Stark  Draper  Laboratory  and 
the  Department  of  Aeronautics  and  Astronautics,  MIT. 

The  air  bearing  table  was  clamped  so  that  the  central  hub  could 
not  rotate.  The  thruster  lines  for  one  of  the  active  appendages  was 
set  to  300  psi  and  the  thruster  was  fired  for  .05  seconds  (50  milli¬ 
seconds).  With  the  appendage  initially  at  rest,  the  firing  of  the 
thruster  was  assumed  to  have  begun  at  time  t  -  0.  Linear  acceleration 
at  the  tip  was  observed  over  the  time  Interval  0  to  5  seconds.  With  a 
sampling  period  of  .005  seconds  (5  milliseconds)  a  total  of  1000 
measurements  were  recorded.  The  data  is  plotted  in  Figure  4.1  below. 
The  soale  factor  for  the  accelerometer  is  5  volts/g  (g  -  32  ft/sec*). 


0  2  4 

TIME  IN  SECONDS 


Figure  4 . 1 

The  noticeably  higher  frequency  (■  14  Hz)  component  of  the  data 
is  a  torsional  mode  of  the  arm  excited  by  the  motion  of  the  thruster 
valve  mechanisms  and  inertial  and  elastic  forces  applied  to  the  tip  of 
the  arm  by  the  nitrogen  supply  hose.  The  opening  or  closing  of  the 
solenoidal  valve  in  the  thruster  generates  an  inertial  force  which 
acts  as  a  torque  on  the  tip  of  the  arm.  Consequently,  torsional  modes 
are  excited.  Also,  in  addition  to  modifying  transverse  bending 
characteristics,  since  the  hose  is  attached  to  the  top  of  the  arm,  its 
horizontal  motion  will  tend  to  generate  torques  which  have  a 
"twisting"  effeot.  Although  the  aocelerometer  is  mounted  at  the 


center  of  the  arm  (and  therefore  on  a  nodal  line  of  the  longitudinal 
torsional  modes,  if  we  assume  vertical  symmetry),  as  the  arm  twists, 
the  accelerometer  picks  up  a  component  of  the  earth's  gravitational 
force.  Since  the  first  torsional  mode  has  a  much  higher  frequency 
than  either  of  the  first  two  flexible  modes  (.75  Hz  and  7.5  Hz,  as 
identified  from  an  FFT  of  the  data)  and  since  it  is  rapidly  damped,  we 
neglected  its  contribution  to  the  accelerometer  signal,  treating  it  as 
white  noise,  and  left  it  unmodeled.  A  detailed  discussion  of  the 
causes  of  the  excitation  of  the  torsional  modes  and  its  effect  on  the 
transverse  bending  characteristics  of  the  active  appendages  can  be 
found  in  [6] . 

The  physical  characteristics  of  the  structure  are  as  follows. 

The  arm  is  made  of  aluminum  and  is  4  feet  in  length,  6  inches  in  width 
and  .125  inches  in  height.  From  this  we  obtain  a.  «  4.0  ft,  p  -  .027 
slug/ft  and  I  -  4.71  x  10-8  (ft)4.  The  theoretically  predicted  value 

Q  O 

for  E  is  15.84  x  10  lb/ (ft)  .  The  mass  of  the  thruster  assembly  was 
determined  to  be  mT  -  .149  slug.  From  the  calibration  table  in  [6], 
we  find  that  a  hose  pressure  of  300  psi  is  equivalent  to  a  force  of 
.297  lb.  We  set  therefore 

,  .  r  0.297  lb  0  S  t  S  0.05 

f(t)  -  { 

0.0  0.05  <  t  £  5.0 

To  serve  as  a  basis  for  comparison,  we  neglected  the  hose  effects 
and  structural  damping  (i.e.  we  chose  "  °H  “  kH  “  0)  and  used 

the  standard  Euler-Bernoulli  model  with  the  parameters  p,  E,  I  and  mT 
and  input  f  as  specified  above  to  generate  the  plot  of  linear 
acceleration  at  the  tip  given  in  Figure  4.2. 
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Figure  4.2 

The  plot  was  obtained  by  integrating  the  initial  value  problem  (3.13), 
(3.14)  with  N  -  4  and  then  using  (4.1)  to  oompute  the  acceleration  at 

3*u  8*uN 

the  free  end.  The  residuals  ( — -(t.Jt) - ( t . JL) )  over  the  time 

8t*  8t* 

interval  [0,8]  are  plotted  in  Figure  4.3.  The  sum  of  the  squares  of 
the  residuals  (at  Intervals  of  .1  seoonds)  was  found  to  be  3.03. 
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Figure  4.3 

Using  the  data  on  the  interval  3.0  to  5.0  (where  the  contribution 
from  the  torsional  modes  has  been  significantly  damped)  with  a 
sampling  period  of  .1  seconds  we  used  our  soheme  with  N  -  4  to  obtain 
optimal  estimates  for  the  coefficient  of  viscosity  oD  and  the  hose 
parameters  mH,  cH  and  kg.  In  the  set  of  runs  we  are  about  to  describe 
the  values  of  E  and  m^  were  held  fixed  at  their  theoretically 
predicted  values.  A  rough  calculation  based  upon  "matching"  the  first 
two  observed  natural  frequencies  of  the  data  with  the  first  two  modal 
frequencies  of  the  model  was  used  to  obtain  a  orude  initial  estimate 


for  the  ratio  kg/mg.  Then,  using  our  scheme  to  minimize  over  the 
parameters  mH  and  kg  only,  we  obtained  the  optimal  values  shown  in 
Table  4.2  below.  Integrating  the  system  (3.13),  (3.14)  over  the  time 
interval  [0,5]  with  m^  and  kg  set  to  the  values  in  the  table  and  « 
Cg  -  0  the  sum  of  the  squares  of  the  residuals  (at  intervals  of  .1 
seconds)  was  found  to  be  .73. 


mH(slug) 

kH(lb/ft) 

.039269 

.339935 

Table  4.2 


Next,  holding  m„  and  kw  fixed  at  the  values  shown  in  Table  4.2,  a 


search  on  cH  was  carried  out  (the  initial  estimate  for  cH  was  taken  to 


be  zero  and  cD  was  held  fixed  at  zero).  Then  using  the  resulting 
values  of  mg,  cH  and  kg  as  initial  estimates,  a  fit  over  all  three 
parameters  was  performed.  The  result  is  shown  in  Table  4.3.  The  sum 
of  the  squares  of  the  residuals  was  found  to  be  .728. 


mH(slug) 

Cgdb-seo/ft) 

kgdb/ft) 

.043431 

.004056 

.351385 

Table  4.3 


Continuing  to  use  the  same  procedure  to  generate  "start  up"  values,  we 


eventually  used  our  scheme  to  search  over  all  four  parameters  cD>  mH 
cH  and  kg  simultaneously  obtaining  the  values  given  in  Table  4.4  and 
the  fit  plotted  in  Figure  4.4.  The  residuals  are  plotted  in  Figure 
4.5.  The  sum  of  their  squares  was  computed  to  be  .70. 


cD(lb-sec/(ft)2) 

mHCslug) 

CgClb- sec/ft) 

kgClb/ft) 

127.40 

.0801 

.007804 

.412977 

Table  4.4 


Figure  4.4 
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Figure  4 . 5 

In  designing  a  controller  for  the  RPL  experiment,  Floyd  in  [6] 
used  model  adjustment  to  tune  a  simple,  undamped,  cantilevered  beam 
with  tip  mass  model  for  the  active  arms  (l.e.  the  arms  with  the  hoses) 
of  the  structure.  He  used  the  following  procedure.  The  air  bearing 
table  was  locked  in  a  stationary  position.  With  the  hose  depressur¬ 
ized,  an  impulsive  force  was  applied  to  the  beam  and  linear  accelera¬ 
tion  at  the  tip  was  measured  and  reoorded.  Based  upon  the  physical 
assumption  that  with  the  hose  depressurized,  the  presence  of  the  hose 
serves  only  to  add  mass  to  the  tip  of  the  arm,  the  parameter  mT  was 


adjusted  so  that  the  first  mode  or  frequency  of  the  model  agreed  with 
the  first  observed  cantilever  mode  (obtained  via  an  FFT)  of  the  data. 
Then,  with  the  hose  pressurized,  the  same  experimental  procedure  was 
carried  out.  This  time  however,  the  modulus  of  elasticity  E  of  the 
beam  was  adjusted  to  compensate  for  the  variation  in  stiffness  which 
results  from  the  presence  of  the  hose.  The  adjusted  values  of  the  tip 
mass,  mT,  and  modulus  of  elasticity,  E,  obtained  by  Floyd  are  given  in 
Table  4 . 5  below . 


mT  (slug) 

E  (lb/(ft)2) 

.254 

17.31  x  108 

Table  4.5 


We  integrated  the  system  (3.13),  (3.14)  using  the  adjusted  values  of 
mT  and  E  given  in  the  table  (and  °d  “  -  cH  -  kH  -  0)  and  obtained 

the  plot  shown  in  Figure  4.6.  The  corresponding  residuals  are 
plotted  in  Figure  4.7.  The  sum  of  the  squares  of  the  residuals  was 
computed  to  be  5.1. 

Starting  with  the  same  basic  model ,  we  used  our  scheme  to 
determine  the  values  of  mT  and  E  which  minimize  the  sum  of  the  squares 
of  the  residuals  over  the  time  interval  [3.0,  5.0]  with  a  sampling 
period  of  . 1  seconds .  Taking  the  theoretically  predicted  values  of  mT 
and  E  (mT  -  .149  slug,  E  -  15.84  x  106  lb/ (ft)2)  as  start  up  values 
for  the  optimization  routine  yielded  the  results  given  in  Table  4.6. 


a 


The  corresponding  fit  and  residuals  are  plotted  in  Figures  4.8  and  4.9 
respectively  below.  The  sum  of  the  squares  of  the  residuals  (over  the 
Interval  [0,53)  was  computed  to  be  .73. 


mT  (slug) 

E  (lb/(ft)2) 

.185 

21.95  x  108 

Table  4.6 
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Figure  4.9 

In  summary,  we  have  seen  that  analysis  of  the  RPL  experimental 
data  can  be  carried  out  in  several  ways  with  a  number  of  different 
models.  Our  techniques  can  be  used  to  provide  reasonable  fits  of  the 
data  to  models  with  or  without  hose  and/or  beam  damping.  Even  if  one 
attempts  to  leave  the  physics  of  the  hose  -  beam  dynamic  interaction 
unmodeled  and  perform  "model  adjustment"  (by  adjusting  the  values  of 
the  tip  mass  mT  and  beam  modulus  of  elasticity  E) ,  our  estimation 
techniques  provide  a  much  better  fit  than  that  obtained  using  "modal 
matohing"  methods  oommon  in  engineering  practioe. 

One  of  the  primary  objectives  of  our  effort  here  was  to 
demonstrate  the  effloaoy  of  our  soheme  and  in  particular,  to  assess 
its  effectiveness  when  provided  with  actual  experimental  data.  While 
we  are  pleased  with  the  results  obtained  for  the  RPL  data,  we  are 
careful  to  point  out  that  to  provide  a  fair  and  oomplete  evaluation  of 


the  usefulness  of  our  models  for  the  RPL  experimental  structure,  a 
more  complete  and  in-depth  study  involving  extensive  experimental  work 
and  statistical  analysis  would  necessarily  be  required. 
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